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Keeping the exact general relativistic treatment of light bending as a reference, we compare the 
accuracy of commonly used approximate lens equations. We conclude that the best approximate lens 
equation is the Ohanian lens equation, for which we present a new expression in terms of distances 
between observer, lens and source planes. We also examine a realistic gravitational lensing case, 
showing that the precision of the Ohanian lens equation might be required for a reliable treatment 
of gravitational lensing and a correct extraction of the full information about gravitational physics. 
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I. INTRODUCTION 

Gravitational lensing is a well-established research sub- 
ject treating the bending of light trajectories by gravita- 
tional fields. Its methodology is traditionally developed 
within the weak field approximation of General Rela- 
tivity, which describes photon trajectories through the 
geodesies equation. In this context, the true position of 
a source in the sky and its apparent position after de- 
flection of the light by a massive body are related by 
the so-called lens equation [l|, |j| (for textbook reviews, 
see [3, 0; Q ) • This relation and its mathematical proper- 
ties have been extensively studied for generic lens models 
and represents the basis of the whole gravitational lens- 
ing theory. 

The lens equation is typically established in the small 
angle approximation, in conjunction with the weak de- 
flection hypothesis. In more recent years, there has been 
a renewed interest in astrophysical situations in which 
the deflection angle is not small. For example, if the 
bending of light emitted by sources near a black hole is 
considered, the deflection angle may reach arbitrary large 
values [g] . For such gravitational lenses the old small an- 
gles lens equation must be obviously revised. Several pro- 
posals for generalized lens equations have then appeared 
in the literature and have been applied to speciflc cases 
0, S, i, 0, UH, E El- Such lens equations principally 
differ among each other for what concerns the variables 
in which they are expressed. However they also lie at 
different approximation levels below the full general rel- 
ativistic description of the photon motion. 

In this work, we review all these lens equations for 
spherically symmetric bodies appeared in the literature 
and introduce a new lens equation that fills a gap in the 
present taxonomy. We also present a detailed discussion 
of the order of magnitude of the errors committed in the 
use of different lens equations, giving a complete inter- 
pretation for their origin. It is worth mentioning that a 
critical review of the approximations leading to the weak 
deflection lens equation has recently appeared [lj|. The 
present work shares the same spirit by testing different 



approximations leading to large deflection lens equations, 
instead. 

The plan of the paper is as follows. In Section II we 
establish the notations by describing the basic geomet- 
ric configuration for gravitational lensing. In Section III, 
we review the exact lens equation by Frittelli, Kling and 
Newman, derived in a fully general relativistic context. 
In Section IV, we discuss the asymptotic approximation 
and all lens equations making use of this approximation. 
We also introduce our new proposal for an improved lens 
equation. In Section V we present a numerical example 
about a realistic gravitational lensing situation, in order 
to compare the different approximations in the lens equa- 
tions previously discussed. Section VI contains a discus- 
sion about the precision needed in gravitational lensing 
observations and the conclusions of the work. 



II. BASIC LENSING GEOMETRY 

Let us consider a spherically symmetric spacetime, 
whose metric is 



ds-" 



A{r)dt^ - Bir)dr^ - C{rf{dd^ + sind^d(f>^). (1) 



Let us put a static source at radial coordinate d^s and 
a static observer at radial coordinate doL- The orien- 
tation of the polar coordinates is chosen in such a way 
that both the source and the observer lie on the equa- 
torial plane d = 7r/2. As a consequence of the spherical 
symmetry, the whole photon motion takes place on this 
plane. We also assume that the metric is asymptotically 
fiat, so that the coordinates (t, r, i), (/>) become Minkowski 
polar coordinates very far from the center (we neglect any 
effects due to cosmological expansion in our discussion). 

A very useful structure for understanding several lens 
equations proposed in the literature is the Minkowski 
space associated with ([T]). This structure can be intro- 
duced by assuming that the metric depends on one or 
more parameters continuously so that, once we tune these 
parameters to zero, the metric becomes Minkowski and 
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FIG. 1: Generic gravitational lensing configuration. Note that only the angles 7, 6 and 9 are well-defined in the curved 
spacetime, whereas all other geometrical quantities need to be defined referring to the Minkowski space associated with the 
original curved space. 



the coordinates (i, r, ??, </>) become Mmkowski polar coor- 
dinates on the whole spacetime. The parameters of the 
metric may be the mass of the lens, the electric charge, 
the scalar charge or any other parameter dictated by any 
kind of gravitational theory. Many definitions need to 
be given by referring to this associated Minkowski space 
obtained by tuning all these parameters to zero. Of 
course, the original curved spacetime and the associated 
Minkowski space coincide in the asymptotically flat re- 
gion. 

A typical geometric configuration of gravitational lens- 
ing is shown in Fig. [I] with a photon emitted by a source 
S", curved by the gravitational field generated by the lens 
L and detected by the observer in O. The observer sees 
the photon at an angle 6 with respect to the optical axis 
OLM, whereas, if there were no lens, he would directly 
observe the source at an angle (3. (3 is thus the first quan- 
tity that is well-defined only in the associated Minkowski 
space, as we must tune the mass of the lens to zero to 
define it. 

The emission direction of the photon is SC , whereas 
the detection direction is CO. In the associated 
Minkowski space, we can also measure the angle between 
these two directions and define it as the deflection angle 
a (we cannot compare directions from different points in 
a curved space). 

It is common practice to define the lens plane and the 
source plane as those planes orthogonal to the optical axis 
passing through the lens and the source, respectively. Of 
course, also these definitions can be given in the associ- 
ated Minkowski space and then extended to the original 
curved space. Then, one can define the distances from 
the observer to the lens plane as Dql, the distance be- 
tween the lens plane and the source plane as D^s a-nd 



the distance between the observer and the source plane 
as Dos- The usual relation 



Dos = D 



OL 



D 



LS 



(2) 



doL 


= DoL 


d-LS 


= D LS / cos^ 


dos 


= Dos/ COS p. 



holds. 

These distances between planes are obviously different 
from the distances between the points O, L and S. We 
have already defined doL and d^s as the radial coordi- 
nates of observer and source in the original curved met- 
ric. Once we report these coordinates in the associated 
Minkowski space, they coincide with the proper distances 
of the observer and the source from the lens. Then, in 
this space, we can establish simple geometrical relations 
with the distances between planes. In particular, we have 



(3) 
(4) 

(5) 



Of course, if the source is very close to the optical axis, 
/3 and 7 are small and the differences between the upper- 
case distances and the lowercase distances is of second 
order in the angles. Therefore, in the classical weak de- 
flection paradigm, the two notions are confused without 
consequences. Here, we must keep them distinct in order 
to avoid confusion. 



III. EXACT LENS EQUATION 

The lens equation is a relation among the source and 
observer coordinates and the angle 6 at which the ob- 
server detects an image of the source S. 



Remaining in a fully general relativistic context, it 
is possible to write down the exact equations govern- 
ing the photon motion and consequently write an exact 
lens equation. This approach has been proposed by Frit- 
telli and Newman in Ref. iTj] and then applied to the 
Schwarzschild lens in Ref. [8| . Later on, it has been gen- 
eralized to all spherically symmetric spacetimes in Ref. 
|la | (see also 16]). 

A photon emitted from a source, deflected by the lens 
and detected by the observer experiences a change in the 
aziniuthal coordinate given by 



<^{J,doL,dLs) = 



doL 



dhS 



ro 



VW)J 



zdr. 



•^Wy^"^ 



(6) 

where J is the specific angular momentum of the photon, 
which is a constant of motion, and tq is the distance of 
closest approach to the lens. These two quantities are 
related by 



J = 






(7) 



The specific angular momentum is also related to the 
angle 9 at which the observer detects the photon. Defin- 
ing the angle by using the scalar product between the 
arrival direction and the direction specified by the opti- 
cal axis OL (see e.g. [l3|)j we find 



— arcsin J\ 



C{doL) 



(8) 



If the observer is very far from the lens, A{doL 
and C{doL) 



dQj^ , so that one recovers the relation 



~ arcsin J/ do L, 



(9) 



which allows to identify J with the impact parameter u 
of the light ray trajectory (except for a speed of light 
factor). In general, however, this identification is only 
approximate (see next section). 

The exact lens equation can finally be written by not- 
ing that the change in the azimuthal coordinate must be 
equal to the difference between the azimuthal coordinates 
of the observer and the source. Replacing all occurrences 
of J and ro in terms of 9 by Eqs. ([7]) and ([S]), we get 

0,i[il 



^i9,d, 



OL,"-LS) 



1- 



(10) 



This equation is the exact general relativistic relation 
between the angle at which the image appears in the ob- 
server's sky and the relative positions of source, observer 
and lens. All other lens equations represent approximate 
forms of this equation. Of course, most lens equations 
are derived under assumptions that are largely satisfied 
in realistic situations. Therefore, their simplifications are 



absolutely welcome, if the induced errors are below ob- 
servational sensitivity or other sources of noise. In the 
following sections we will introduce some lens equations, 
describing their approximations. 



IV. APPROXIMATE LENS EQUATIONS 

The most popular approximation is what can be called 
as asymptotic approximation, which amounts to saying 
that the source and the observer are in the asymptotic 
flat region of the spacetime. Quantitatively, this is ex- 
pressed by requiring doL, d^s ^ fg, where Tg is the grav- 
itational radius of the lens, i.e. the typical scale in the 
curved spacetime metric controlling the range of the grav- 
itational field (it is 2GM /(? in the Schwarzschild metric). 
As a first consequence, Eq. ^ holds, so that the angu- 
lar momentum of the photon J can be identified with 
the impact parameter u of the initial trajectory and the 
impact parameter u of the final trajectory, as depicted 
in Fig. [T] As a second consequence, the deflection angle 
can be calculated as 



a{u) = $( J — u, oo, cx)) ■ 



(11) 



which is the azimuthal shift of a photon incoming from 
infinity with impact parameter J and escaping to infinity. 

Once this approximation is accepted, it is possible to 
establish a relation among 9 and the relative positions of 
source, lens and observer using pure Euclidean geometry. 
In fact, all distances are defined in the asymptotic space- 
time and the only input from General Relativity is the 
precise expression of the deflection angle as a function of 
the impact parameter u. 

Note that the asymptotic approximation does not nec- 
essary imply the small angle approximation. In fact, it 
just assumes that the source and the observer distances 
are much larger than the gravitational radius of the lens. 
This does not prevent the impact parameter to be very 
large. 

In the following subsections, we will introduce several 
approximate lens equations appeared in the literature, 
expressed in terms of different quantities. 



A. Ohanian lens equation 

The first lens equation we introduce is due to Ohanian 
^f\ , who proposed it in a study of gravitational lensing by 
a Schwarzschild black hole for arbitrary defiection angles. 

Thanks to the asymptotic approximation, we can use 
Euclidean geometry to relate the various quantities and 
store all the relativistic input in the angle a. Let us define 
the angle C = OSL. Considering the triangles OLS and 
OSC, we can write down the relations 



P + C + TT-J^TT 

{9^(3) + {9~C)+TT-a 



(12) 
(13) 



Summing up these two equalities, we obtain [T^ (see also 

e + e-a = j. (14) 



The angle 9 can be expressed in terms of 9 recalling 
that the impact parameter of the incoming trajectory u is 
equal to the impact parameter of the outgoing trajectory 
u (in the asymptotic approximation). Therefore 



fdoL 



arcsm — — sm ( 



(15) 



Once Eq. p5|) is used in Eq. [TH we obtain the lens 
equation as a relation involving the detection angle 9, 
the distance between the observer and the lens doL, the 
distance of the lens to the source dis, the source position 
angle 7 and the deflection angle a, which is a function of 
9 through u = doL sm9. 

The original form proposed by Ohanian [9[ was actu- 
ally written by replacing 9 and 9 by their small angle 
approximations u/dpL and u/dLs respectively. This is 
an additional approximation, with respect to the asymp- 
totic approximation, which we will not consider here. 



B. VNC lens equation 

Another very simple expression for the lens equa- 
tion was proposed by Virbhadra, Narasimha and Chitre 
(VNC) in a paper studying the role of the scalar field in 
gravitational lensing [lOJ. It is just the result of the sine 
theorem applied to the triangle OCS 



sm{9 — P) = - — sin a, 
dos 



(16) 



where dcs is the distance between the source S and the 
intersection point C between the incoming and outgoing 
trajectories. Eq. pB]) is valid for arbitrary deflections 
and impact parameters. As for the Ohanian lens equa- 
tion, its only approximation with respect to the exact lens 
equation pop is the asymptotic approximation, which is 
necessary to express the deflection angle a. 

However, the lens equation P^ is not a closed relation 
among the detection angle 9 and the relative positions of 
observer, lens and source. In fact, it is expressed in terms 
of the position of point C, which is unknown a priori and 
needs to be estimated in some way. One possibility is 
to approximate dcs by dis, which is reasonable if the 
source is at distances much larger than the impact pa- 
rameter. However, this additional approximation would 
spoil the effectiveness of the lens equation in its original 
form. Therefore, the VNC lens equation is difficult to 
use in practical applications, since it demands an inde- 
pendent knowledge or at least an estimate of the position 
of the point C. 



C. Virbhadra and Ellis lens equation 

With their outbreaking paper about the possibility of 
observing higher order images around the black hole at 
the center of our Galaxy [ll| , Virbhadra and Ellis have 
attracted a great attention on gravitational lensing be- 
yond the weak deflection approximation, inspiring new 
vitality in black hole gravitational lensing. They have 
also proposed a new lens equation that has become very 
popular in the scientific literature. Their equation is writ- 
ten in terms of the distances between source, lens and 
observer planes. 

Starting from the relation among the segments 



MS = MI - SI, 
we can write the relation 



(17) 



Dos tan /? — Dos tan 9 — Dcs [tan 9 + tan(Q; — 9)] , 

(18) 
where Dcs is the distance between point C and the 
source plane (namely the length of the segment CH). 
Even in this case, the lens equation is expressed in terms 
of the position of point C, which should be estimated in 
some way. The proposal by Virbhadra and Ellis is to as- 
sume that C lies on the lens plane, so that Dcs — Di^g. 
Then the final form of the lens equation is 

Dos tan /3 = Dos tan 6 — Dls [tan 9 + tan(a — 9)] , 

(19) 
At the cost of an additional approximation, the lens equa- 
tion is put in a form very easy to use in many astrophys- 
ical applications, as it is finally expressed in terms of the 
positions of lens, source and observer. We will discuss 
the error introduced by this approximation in Section FVl 
while in Section [IV Fl we will present an improvement of 
this equation which avoids this approximation and makes 
it equivalent to the Ohanian lens equation. 



D. Dabrowski and Schunck lens equation 



In a paper studying gravitational lensing by boson 
stars [i2], Dabrwoski and Schunck realized the difficul- 
ties of using the VNC lens equation ITBl) and derived the 
alternative lens equation 



1(0-/5) 



dLS 

dos 
[tan 9 - 



■ cos f cos < arcsm 



tan(a - 9)] . 



dos 

dLS 



sin/3 



(20) 



This can be obtained from the Virbhadra and Ellis 
lens equation p9)) replacing the distances between planes 
Dos and D^g by the distances between objects using 
Eqs. ^ and (O and noting that 



dLssin-f = dos sin (3 
in the associated Minkowski space. 



(21) 



In practice, the equation ((20|) by Dabrwoski and 
Schunck is the same as Eq. p^ by Virbhadra and El- 
hs but expressed in terms of distances between objects 
instead of distances between planes. Of course, hke the 
former equation, it contains the additional approxima- 
tion that C lies on the lens plane. 



E. Bozza and Sereno lens equation 

The Ohanian lens equation P^ has the advantage of 
being very simple and being the closest relative of the 
exact lens equation, since it only contains the asymp- 
totic approximation and makes no additional assump- 
tions. However, in several astrophysical applications one 
may prefer to have a lens equation directly written in 
terms of the angle /3 rather than 7. In fact, (3 is an angle 
with vertex in the observer and thus directly connected 
to coordinates in the observer's sky, whereas 7 is an angle 
with vertex in the lens, thus being less close to observ- 
ables. The relation between the two angles is given by 
Eq. (HI]). Therefore, if we take the sine of equation (fH)) . 
using Eq. p^ and re-ordering terms we get 



(io5 sin /3= do^ sin6'cos(a — 



dls 



"-OL '• 



d), (22) 



which first appeared in Ref. |13| . Though being more 
complicated than Eq. (fT4|) . this lens equation has the 
advantage of being directly expressed in terms of the an- 
gle (3 and thus preferable for applications in which this 
angle is directly involved. Otherwise, it is completely 
equivalent to the Ohanian lens equation. Note that the 
distances involved in this equation are distances between 
the objects and not between their geometrical planes, as 
in the lens equations of Sections IIVCI and IIVFI 

One may note that this lens equation is expressed in 
terms of the same variables appearing in the equation by 
Dabrwoski and Schunck (P0|) . However the two equations 
look different. Apart from the ordering of the terms, Eq. 
([2OI1 contains the additional approximation that C lies 
on the lens plane, whereas Eq. ([221) does not. Therefore, 
Eq. ((22I) represents an improved version of Eq. ((20)) . 
The same difference will arise in the new equation to 
be presented in the next section and the equation by 
Virbhadra and Ellis. 



Starting again from the Ohanian lens equation (|14p . 
we can solve it in terms of 9 and plug it in Eq. (fT5|) . 
Recalling Eq. ^ and the relation between 7 and (3 (pij) . 
we can put the lens equation in the form 



Dos tan f3 — 



Dql sin6> - Dls sm{a 
cos{a — 9) 



(23) 



which represents the improved version of Eq. (|19p by 
Virbhadra and Ellis, in the same sense as Eq. ([22)) by 
Bozza and Sereno represents the improved version of Eq. 
([^0)) by Dabrowski and Schunck. In Section V we make a 
thorough discussion and comparison of the two equations. 



G. Small angles lens equation 

We finally recall the classical lens equation obtained in 
the hypothesis of small angles a,9,f3<^ 1. This hypoth- 
esis has nothing to do with the weak deflection approx- 
imation, since the deflection angle is typically expressed 
in powers of rg/u, whereas the corrections to the small 
angles approximation are expressed in powers of 9, a, (3. 
So it makes sense to consider an exact deflection angle 
while performing the small angles approximation in all 
trigonometric functions. 

The small angles lens equation can be obtained from 
any of the equations j TO t . ([20)) . ([22)) or ^. The well- 
known result [U, [3, S m, m is 



(3 = 9- — — a. 
Dos 



(24) 



The small angles approximation is the most rude ap- 
proximation that can be done on the lens equation. Nev- 
ertheless, it is a useful reference approximate equation 
lying on the other extremum of the approximation lad- 
der with respect to the exact lens equation ([TU)) . with all 
other lens equations staying in the middle steps between 
the two equations. 

It is interesting to note that a similar equation can be 
deduced when the deflection angle is very close to mul- 
tiples of 27r. This case corresponds to higher order im- 
ages generated by photons performing one or more loops 
around the lens before emerging. The only difference is 
that the deflection angle a in Eq. ([M)) must be replaced 
by a — 2ri7r, where n is the number of loops. 



F. A new improved lens equations between planes 



H. Conclusions on approximate lens equations 



Eq. (|19p represents a very useful lens equation ex- 
pressed in terms of distances between the observer, lens 
and source planes. However, as we have explained in Sec- 
tion [IVCl it is derived under the additional assumption 
that the intersection point C between the incoming and 
outgoing ray trajectories lies on the lens plane. It would 
be desirable to have a lens equation expressed in terms of 
the same quantities without this additional assumption. 



To summarize this section, we can say that the asymp- 
totic approximation allows the use of the relations of 
Euclidean geometry, confining General Relativity to the 
derivation of the precise expression of the deflection angle 
in terms of the impact parameter. 

In order to write the lens equation in terms of 9, a 
and the relative positions of observer, source and lens, 
an additional relation is needed between the incoming 



and the outgoing branches of the photon trajectory. The 
lens equation by Ohanian p^ . and the related equations 
([22| by Bozza and Sereno and ([23|) presented in this pa- 
per, use the equality of the impact parameters of the in- 
coming and outgoing trajectories, expressed by Eq. (fT5)) . 
This equality descends from the time-reversal symmetry 
of the photon geodesic in General Relativity along with 
the asymptotic approximation. 

On the other hand, the lens equations by Virbhadra 
and Ellis (fT9|) and Dabrowski and Schunck (|20|) use an 
alternative relation between the incoming and outgoing 
trajectories by imposing that the intersection point C lies 
on the lens plane. In the next section we will quantify 
the accuracy of these lens equations in a realistic astro- 
physical situation. 



V. LENS EQUATIONS AT WORK IN A 
REALISTIC EXAMPLE 

As explained in the previous section, in the treatment 
of gravitational lensing by spherically symmetric lenses 
with arbitrary deflection angles, we can distinguish an 
exact lens equation (jlOp , a family of lens equations with 
the asymptotic approximation only p^ . ^F^ and (P5)l . 
and a family of lens equations with the additional geo- 
metrical approximation that the intersection point C lies 
on the lens plane (fT9|) and ((20|) . For simplicity, we shall 
refer to these two families of approximate lens equations 
as the Ohanian family and the Virbhadra and Ellis fam- 
ily, respectively. Furthermore, we have the small angles 
lens equation, in which all trigonometric functions are 
approximated by their first order expansions. Finally, 
we can also consider the effect of the weak deflection ap- 
proximation on the deflection angle, by retaining only the 
lowest order term in the weak deflection expansion. We 
will not consider here the VNC lens equation (fTB|) since 
it is not expressed in terms of the relative configuration 
of source, lens and observer but demands additional in- 
formation on the position of point C . 

As pointed out in Ref. [ll|, the best candidate for 
the observation of gravitational lensing in the regime of 
large deflections is the black hole at the center of our 
Galaxy, named Sgr A*. Its mass is estimated to M — 
3.6 X 10^ Mq and its distance is Dql = 8 kpc [11]. The 
Schwarzschild radius is therefore 



2GM 



lO^^m. 



(25) 



In order to test the accuracy of the different lens equa- 
tions, we consider the source plane to be at distance 
Dls — lOOOrg from the black hole. At this distance, 
the source is just at the margin of the accretion disk of 
the black hole, so it may be a star spiralling towards the 
black hole (the star S2 is known to reach a distance of 
1700rg at its periapse [1^, |23|) or a bright spot of hot 
material rotating around the black hole. Apart from the 
physical meaning of a source at such distance, this choice 



has been made with the aim of putting in better numer- 
ical evidence the differences in the performances of the 
lens equations, since the expected order of magnitude of 
the errors scales with powers of Vg/ D^s and Vg/ Dql- 

We are going to study the position of the images as 
a function of the position of the source, keeping the dis- 
tance to the source plane Dqs fixed. Therefore, we will 
compare the results of the exact lens equation pO|) with 
the results of the Virbhadra and Ellis (VE) lens equation 
(|19p . the improved lens equation (|23p . which is equivalent 
to Ohanian lens equation (wc shall refer to it as the OB 
lens equation), and the small angles lens equation (|24p . 
The last three equations are already expressed in terms 
of /3, Dls and Dqs, whereas the exact lens equation p^ 
is expressed in terms of 7 and d^s- We will use Eqs. (^J) 
and ((4]) to obtain these quantities as functions of /3 and 
Dls, referring to the associated Minkowski spacetime. 

The azimuthal shift $ is calculated assuming that the 
spacetime metric around the black hole is Schwarzschild. 
Therefore we have 

A(0 = (l-^); S(r)=(l-^)"'; C{T)^r^. 

(26) 
We remind that the detection angle 6 is related to the 
angular momentum of the photon J by Eq. ([S]) . In the 
approximate lens equations (fT9| . (|23p and ((24l) the de- 
flection angle is obtained by Eq. (fTTj) as a function of 6 
through Eq. ^ . For any given value of the source posi- 
tion j3, we find the corresponding position of the images 
9 numerically in all lens equations (see Section IV CI for 
details). 

A Schwarzschild black hole generates two infinite se- 
quences of images on each side of the source '6|. The 
outermost pair of images is made of the classical primary 
and secondary images, which can also be described in the 
classical weak deflection limit. Section IV Al ls devoted to 
their study. The second pair of images is generated by 
photons performing one complete loop around the black 
hole before reaching the observer. These higher order im- 
ages are typically very faint, but their detection could be 
very important for a confirmation of General Relativity. 
We examine them in Section fV Bl 



A. Primary and secondary images 

For any values of /3, we look for solutions of the lens 
equation such that 7r/2 < $(0, doL,dLs) < 3n/2. In this 
interval we always have one solution for any values of (3. 

In Fig. [2] we show the exact position of the image for 
f3 in the range [—10 mas, 10 mas]. For positive /3 the im- 
age is on the same side of the source (primary image), 
whereas for negative f3 the image is on the opposite side 
(secondary image). We see that for large values of /3 
the position of the primary image tends to coincide with 
the source position itself {9 ~ /?). The secondary im- 
age, instead, becomes closer and closer to the black hole. 



4^ 




FIG. 2: The exact position of the image as a function of the 
source position /3 in a linear-log plot. The dashed line is the 
shadow of the black hole 0,nin ■ For /3 > the plot represents 
the primary image, for /3 < it represents the secondary 
image. 



The minimum angle 9min represents the border of the so- 
called shadow of the black hole. This angle is obtained 
by the minimum angular momentum through Eq. ([5]) 
with Jmin = 3\/3rg/2. For the black hole in Sgr A*, we 
have Ojnin = 23 /ias. The value of 9 for (3 = represents 
the radius of the Einstein ring, which for our geometrical 
configuration is 9e = 404 /ias. 

Fig. [21 has been obtained with the exact lens equation. 
However, as anticipated before, we expect errors at most 
of the order rg/ D^s using the approximate equations. 
So, the difference cannot be appreciated by superposing 
the solutions obtained with different lens equations. We 
choose, instead, to plot the relative error in the image 
positions with respect to the exact lens equation. These 
are defined as 



All lens equations have the same accuracy for very 
large (3, when gravitational lensing becomes negligible. 
However, as (3 approaches 9e the error of all equations 
increases, save for the OB lens equation. Noting that the 
plot is in a logarithmic scale, we see that 6ve tends to be 
of the order 10^"^, whereas Sqb drops to 10^^, proving 
to be much more accurate. Going to /? < 0, we see that 
the error stays nearly constant for the VE lens equation, 
whereas it continues to decrease for the OB lens equation 
to very tiny values. Note the change in sign in 6ve for 
positive /?, which can be deduced from the negative spike 
in the logarithmic plot as 5ve crosses zero. This signals 
the fact that the sources of error dominating for large (3 
and small /3 are different and have opposite signs. Sur- 
prisingly, the small angles lens equation SA has an error 
comparable to the VE lens equation in the most inter- 
esting lensing region (|/?| < 9e), corresponding to a very 
good alignment of the source with the lens. The VE lens 
equation performs better than the SA equation at large 
negative f3, where the small angles approximation is no 
longer tenable. On the other hand, the weak deflection 
approximation WD always stays at larger errors. 

Now let us find the interpretation of the errors of the 
considered lens equations. We start from the OB lens 
equation, which only contains the asymptotic approxi- 
mation. This lens equation is completely equivalent to 
the Ohanian lens equation, so we can use this equation 
for the evaluation of the error and translate the results 
in terms of the OB equation. 

If we expand the azimuthal shift $ ([6]) in inverse powers 
of DoL and D^s, at zero order we obtain the deflection 
angle a + tt. After this, we obtain some terms that re- 
produce the power expansions of the geometrical terms 
—9 and —9, which are explicitly present in the Ohanian 
lens equation (fT4|) and finally, the first term that is not 
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where the subscript "ex" refers to the result of the exact 
lens equation (|10p . "OB" refers to the improved equation 
presented in this paper (P5)) derived from that by Oha- 
nian, "VE" refers to the equation by Virbhadra and Ellis 
P^ . "SA" refers to the small angles equation (|M|) . and 
"WD" refers to the small angles lens equation with the 
Einstein approximation a ~ 2rgl Dol9 for weak deflec- 
tion angles. Apart from 9wd and 9ex^ in all remaining 
cases the position of the images is computed picking a 
from Eq. (|lip and relating u to 6* by Eq. ([8]), as men- 
tioned before. 

Fig. [3] shows a comparison between the accuracy of 
the four equations. 
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FIG. 3: Relative error in the position of the images for the 
weak deflection lens equation {WD), the small angles lens 
equation {SA), the Virbhadra and Ellis lens equation {VE) 
and the improved Ohanian lens equation {OB) in a log scale. 
The dashed vertical lines bound the region with |/3| < 6e, in 
which gravitational lensing is mostly interesting. 
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FIG. 4: Comparison among the relative error in the posi- 
tion of the images for the Virbhadra and Ellis lens equation 
(VE) and what we expect from the asymptotic approximation 
(asym) and the approximation that the intersection point C 
lies on the lens plane (cross). 



present in the Ohanian lens equation is 
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Taking this term to the right hand side of the Ohanian 
lens equation, we can interpret it as an effective change 
in 7. Then we can estimate the corresponding change in 
simply by 






(32) 



Evaluating the derivative of 6 numerically, and taking 
for Ajasym the expression in Eq. pip , we exactly re- 
produce the error eos in the OB lens equation, as can 
be appreciated from Fig. ID In this figure we have also 
shown the VE error, already presented in Fig. [3l together 
with the error deriving from the approximation that C 
lies on the lens plane. This error can be simply calcu- 
lated by taking the differences of the right hand sides of 
Eqs. (PD and ^ 
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(33) 
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We can explicitly see that the error due to the approx- 
imation that C lies on the lens plane dominates the error 
due to the asymptotic approximation in the calculation 
of the secondary image (/3 < 0) and for the primary im- 
age up to moderately large (3. For P ^ 9e, the error 
due to the asymptotic approximation becomes more rel- 
evant than the error due to the position of the intersec- 
tion point. However, as already mentioned before, the 
last regime is the least important for gravitational lens- 
ing. 



The error induced by the small angles approximation 
can be obtained by considering the next-to-leading order 
in the power expansions of the trigonometric functions 
present in Eq. (P5)) . It amounts to 
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To conclude this subsection, it is interesting to give 
approximate estimates of the order of magnitude of the 
errors in terms of the geometric distances involved in the 
problem in different limits. For example, we can consider 
the two limits P :s> 9e (primary image) and — /3 :^ 9e 
(secondary image). In the first case we have ~ /3, in the 
second case 9 ~ —2rgD ls / {DolDqsP) from the weak 
deflection approximation. Taking also Dqs — Dql^ we 
have the simple estimates 
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which reproduce the curves in Figs. [3] and |4] fairly well 
outside the lensing zone and are useful to realize the order 
of magnitude of the errors. In particular, we can appreci- 
ate that the error due to the position of the intersection 
point has opposite sign with respect to the others. 

Finally, a very interesting limit is the case P — 0, cor- 
responding to a source perfectly aligned with the black 
hole and the two images merged into an Einstein ring of 
radius 9. The errors in the estimate of the radius of the 
Einstein ring are 
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(42) 
(43) 
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where the error due to the weak deflection approximation 
comes from the first term neglected in the weak deflection 
expansion (see for example Refs. |2ll. 123. |23| ) . 



These relations are very useful to understand the accu- 
racy of the various lens equations in the most interesting 
regime, that is when the gravitational lensing images are 
most prominent and eventually form an Einstein ring. 
The weak deflection lens equation has an error of the 
order (rglDhs)^- In the case under examination, this 
would translate into an error of 6.6 [is& in the radius of 
the Einstein ring, which is Be — 404 /^as in the situa- 
tion imagined in this calculation. The small angles lens 
equation employed with the exact deflection angle has 
an error of order (r^/Dig), which would amount to 0.07 
/xas. The VE lens equation has an error of the same order 
of magnitude, but with a slightly larger numerical coef- 
ficient, thus leading to an error of 0.2 /ias. Finally, the 
OB lens equation has an error of order [vg/DLsY , which 
/zas. 
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FIG. 5: Fractional displacement of the first pair of higher 
order images from the shadow border as a function of /3 as 
defined in Eq. (|46)>. 



B. Higher order images 



An analysis similar to Section IV Al can be repeated 
for the higher order images, whose position and mag- 
nification has been analyzed in Ref. [ll| for the case 
of the black hole in the center of our Galaxy, and then 
revisited in many papers with different methods. The 
first pair of higher order images can be found in the in- 
terval 57r/2 < ^{9,doL,dLs) < 77r/2, corresponding to 
photons performing one loop before reaching an observer 
on the other side of the lens with respect to the source. 
Since these images appear very close to the shadow bor- 
der 9mim we think it is more instructive to discuss the 
fractional displacement of the image from the shadow 
border, defined as 



(46) 



This quantity is shown in Fig. [5] for the geometric con- 
figuration examined in this paper. The plot has been 
calculated using the exact lens equation (fTU)) . For /3 > 
the curve represents the fractional displacement of the 
positive parity image on the same side of the source and 
for /3 < the curve represents the fractional displace- 
ment of the negative parity image on the opposite side. 
For f3 — 0, we obtain the displacement of the first order 
Einstein ring 9e.i from the shadow border. 

In order to compare the performances of the different 
lens equations, in Fig. [S] we plot the relative error in 
this fractional displacement with respect to the exact lens 
equation 
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FIG. 6; Errors in the estimate of the fractional displacement 
of the first pair of higher order images from the shadow bor- 
ders. SA refers to the small angles lens equation, VE to 
the Virbhadra and Ellis lens equation, OB to the improved 
Ohanian lens equation. The dashed lines bound the range in 
which 1/3 1 < 9e,i, where 6e,i is the angular radius of the first 
order Einstein ring. 



Again we see that the OB lens equation has much lower 
errors in the estimate of the position of the images than 
any other lens equations. The VE lens equation performs 
slightly better than the small angles equation for \/3\ > 
9e,i, whereas the opposite occurs in the good alignment 
regime \(3\ < 9e,i- It is interesting to note that there is a 
value for P such that a = 27r and [3 — 9. In this particular 
point, all lens equations are equivalent and are affected 
by the same small error, determined by the asymptotic 
approximation. 

Similarly to the case of the primary and secondary im- 
ages, we can make analytical estimates of the errors in 
the positions of the higher order images. We follow the 
same procedure described in the previous subsection and 
use the analytical approximation 



The SA lens equation, in this case, is Eq. ([M|l with a 
replaced by a — 27r, as explained in Section [I VGI 



e ~ esDL = 216(7 - 4^3)6-^-+^, 



(50) 
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derived in the strong deflection limit [g, M, ^M i ^^ order 
to have fully analytical results. Then, the errors in the 
displacement of the first order Einstein ring are 
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(51) 
(52) 
(53) 



Again the error in the VE lens equation is of the same 
order as the error in the small angles lens equation with 
a slightly larger numerical coefficient. 



C. Notes on numerical integration 

In the previous subsections we have compared the im- 
ages obtained by solving different lens equations. The rel- 
ative error has proved to be quite small, particularly for 
the OB lens equation, whose error goes down to 10~^^ in 
some plots. In order to investigate such tiny differences, 
we must push the numerical precision of our calculations 
sufficiently far. 

The crucial step in the computations is represented by 
the evaluation of the integral ([6]), which contains an in- 
tegrable singularity at tq. We have used the NIntegrate 
routine by Mathematica with the default method (adap- 
tive Gaussian quadratures with error estimation based 
on Kronrod points [25|). In order to reduce the errors, 
we have worked with the variable z, defined by 



ra 



1-2' 



(54) 



z ranges from to 1 as r ranges from tq to 
Schwarzschild metric, the integrand in Eq. ([( 
in terms of z becomes 



-cx). In the 
|) expressed 



y/zy/2ro - 3 + (3 - ro)z - z^ 



(55) 



We have checked that increasing the precision goal and 
the working precision makes the numerical results con- 
verge to an asymptotic limit, with a controllable precision 
(we have reached a 10~^^ precision). Furthermore, we 
have also re-written the integral as an ordinary differen- 
tial equation to be solved by the NDSolve routine, which 
switches between a non-stiff Adams method and a stiff 
Gear backward differentiation formula method. Solving 
the differential equation, we get the primitive function, 
which can be evaluated at both ends of the domain. In 
this way, we have double-checked our results obtained 
by NIntegrate finding that the two methods converge to 
the same result. In particular, we find that the differ- 
ential equation method converges more slowly than the 
numerical integration with the Gaussian quadratures, as 



the precision goal is increased. For this reason we tend 
to give a preference to the NIntegrate routine. The re- 
sults presented in the plots are calculated with a 10^^^ 
precision. 

The accuracy in the numerical calculation is finally 
confirmed by the very good agreement with the analyti- 
cal estimates for the errors at leading order. 



VI. DISCUSSION AND CONCLUSIONS 

Of course, an approximation can be valid or not de- 
pending on the precision of the observational instrumen- 
tation, the level of environmental noise and the accuracy 
needed for the extraction of the interesting information. 

The best resolution we can expect to reach in a 
reasonably near future is th at of the MAXIM project 
( http://maxim.gsfc.nasa.gov| , which amounts to a frac- 
tion of fj,as. Then, in the example considered in this sec- 
tion, the weak deflection lens equation would not be able 
to ensure such a precision. The VE and the small angles 
lens equations would be marginally adequate, whereas 
the Ohanian lens equation would be largely sufficient, 
without need to resort to the exact lens equation. So, in 
order to exploit the resolution of future instruments, it 
might be necessary to use the Ohanian lens equation or 
one of its variants. 

Environmental noise can be provided by scattering and 
absorption processes affecting the photons involved in 
gravitational lensing. Such processes should be very ef- 
fective in the environment of the supermassive black hole 
at the Galactic center. However, they act in a statistical 
way and a good average process should be able to restore 
the original information regarding gravitational lensing. 
More subtle are the effects that introduce unwanted sys- 
tematics, such as gravitational lensing by secondary ob- 
jects. However, the gravitational lensing optical depth 
towards the center of the bulge is of the order 10~^ ^51], 
which indicates that lensing or microlensing effects by 
stars or other compact objects on the line of sight is abso- 
lutely negligible (see also [2^1 )• In the absence of relevant 
systematics, the position of the centroids of the images 
can be used for a clean reconstruction of the gravitational 
lensing event at the angular resolution level reached by 
the instruments. 

Finally, we note that there are many theoretical rea- 
sons to require a reconstruction of a gravitational lensing 
event as precise as possible. As anticipated, the study of 
higher order terms in the deflection angle requires a pre- 
cision better than {rg/ DlsY^"^ ^^ second PPN order and 
{rg/ Dls) at third PPN order. Therefore, if one wants 
to compare the results of PPN formalism with different 
black hole metrics, it is mandatory to use a lens equation 
that guarantees the necessary accuracy. In this respect 
the VE lens equation does not represent an improvement 
with respect to the small angles approximation, as it is 
affected by an error of the same order of magnitude. 

In this work we have compared the level of accuracy 
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of several lens equations, keeping the exact general rel- 
ativistic lens equation introduced by Frittelli, Kling and 
Newman [3, Q as the reference equation. We have shown 
that the Ohanian lens equation [9] and its close relatives 
[13l llSj are the best approximations of the exact lens 
equation, in that they adopt the asymptotic approxima- 
tion only. Other lens equations are not expressed in terms 
of relative positions of source, lens and observer [lOJ or 
introduce additional approximations [111 . Il2l | . 

Furthermore, we have presented a new formulation of 
the Ohanian lens equation in terms of the distances be- 
tween the observer, lens and source planes, which fills a 
gap in the lens equation zoo. As shown by a numerical ex- 
ample describing a realistic gravitational lensing event by 
the black hole at the center of our Galaxy, such an equa- 



tion represents a noteworthy improvement with respect 
to previous commonly used lens equations expressed in 
terms of the same quantities. 
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